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Rare decays of B mesons, such as B ^K*Y and B K^*'^£'^£^ are loop suppressed in the Standard 
Model and sensitive to new physics. The final state meson in heavy-light decays at large recoil has 
sizeable momentum in the rest frame of the decaying meson. To reduce the resulting discretization 
errors we formulate the nonrelativistic heavy quark action in a moving frame. We discuss the 
perturbative renormalization of the leading order heavy-light operators in the resulting theory 
which is known as m(oving)NRQCD. 

We also present radiative corrections to the NRQCD action computed using automated lattice 
perturbation theory. By combining this technique with high-beta simulations in the weak coupling 
regime of the theory higher order loop corrections can be calculated very efficiently. 
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1. Introduction 

Exclusive decays of B mesons can be used to study the heavy flavour sector of the Standard 
Model. The magnitude of one of the least well known CKM matrix elements, Vub, can be ex- 
tracted from measurements of the decay B niv. In addition, rare processes such a& B ^ K*y 
and B K^*Hl provide an excellent opportunity to constrain new physics models. Recently, the 
experimental uncertainties in these processes have been reduced to below the 5% level [|l|, |^, |3|]. 
Theoretical predictions of comparable precision are necessary to extract fundamental Standard 
Model parameters and to test new physics models. 

The calculation of hadronic heavy-light form factors involves several challenges: on currently 
available lattices, the Compton wavelength of the heavy quark is smaller than the lattice spacing and 
the relativistic quark action can not be used directly. Instead we work with a nonrelativistic effective 
heavy quark action; discretization errors are under control and can be reduced systematically by 
improving the action. If q = ps — Pf is the momentum transfer between the initial and final state 
mesons then can be very small. For B — > niv most experimental data [Q, ^ |6|] comes from the 
large recoil region, and q^ is zero for the radiative rare decay B ^ K*y. This leads to substantial 
discretization errors in the final state meson due to its large momentum. These can be reduced by 
discretizing the nonrelativistic action in a moving frame of reference. The resulting formalism is 
known as m(oving)NRQCD [0, |8|]. 

The heavy quark action used in this work is correct up to ^(1/m^, v^gj), where Vrei is the 
relative velocity between the two quarks in a heavy-heavy meson. One source of systematic errors 
is the mismatch between the ultraviolet modes in the continuum and lattice mNRQCD. In heavy- 
heavy spectrum calculations radiative corrections are of comparable size to relativistic corrections 
[^]. Typical loop momenta are of the order of the heavy quark mass, where the strong coupling 
constant is small. In this work we calculate radiative corrections to the mNRQCD action and heavy- 
light operators. The calculations are carried out at one loop in mean field improved diagrammatic 
perturbation theory. Due to the complexity of the action the Feynman rules are generated using an 
automated expansion algorithm [ |Tq| , |i|, 12] and phase space integrals are solved numerically with 



the adaptive Monte Carlo integrator Vegas [|T^. 

In addition, fundamental parameters of the Standard Model, such as the heavy quark mass 
in the MS scheme, m^'^, can be extracted from lattice calculations if their perturbative expansion 
is known to high orders both in the continuum and on the lattice. In this work we use a mixed 
strategy for calculating higher order loop corrections to the heavy quark self-energy: by simulating 
the quenched theory in the weak coupling regime (corresponding to large values of the inverse 
couphng j8) and fitting the results to an expansion in as, the gluonic higher order loop corrections 



can be obtained [|14|]. This has to be supplemented with the calculation of the fermionic contribution 
in diagrammatic perturbation theory. At one loop order the fermionic vacuum polarisation does not 
contribute and we compare results from diagrammatic perturbation theory and high-j8 simulations. 

1.1 Effective theories for heavy quarks on the lattice 

After integrating out fluctuations at the b quark scale, the Lagrangian of the nonrelativistic 
heavy quark theory can be written as an expansion in the inverse heavy quark mass. A mNRQCD 
action which is correct to ^(1/ni^, v^gj) is derived in detail in [^. On the lattice it can be written as 
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p p + k p P p 

Figure 1: One loop corrections to the heavy quark self-energy. 



. (1.1) 



The lowest order kinetic term is Hq = —iv ■ — ^A^^' — Ay^^ j / {lym) where v is the frame velocity 

and A^, A(2) and Aj ' are first and second order covariant finite difference operators. 5H contains 
higher order corrections in 1/m and operators which reduce discretization artifacts. 

2. Perturbative renormalization of m(oving)NRQCD 

Renormalization constants of kinetic terms in the mNRQCD action can be obtained from the 
heavy quark self-energy, which at one loop is given by the two diagrams in Fig. We include 

(2) 

mean field corrections by dividing each link in the action by its mean value uq = \ — asu:^ + 

2.1 Renormalization parameters 

We calculate the one loop corrections to the zero point energy Eq, wavefunction renormal- 
ization^ 5Z|^, heavy quark mass 5Z,„, frame velocity dZ^ and the energy shift between lattice 
mNRQCD and continuum QCD 5Cy for a range of frame velocities, see Fig. ^ The calculation is 
performed in an infinite volume with a small gluon mass A as an infrared regulator. All renormal- 
ization parameters are obtained from the self-energy by taking appropriate derivatives with respect 
to the external quark momentum 1 15, |l^]. Of particular interest is the renormalization of the exter- 



nal momentum 5Zp, which in the continuum is protected by reparametrization invariance. On the 
lattice, it is small for not too large frame velocities. The gluon action is Symanzik improved in all 
calculations presented here. After including mean field corrections, the magnitude of the one loop 
corrections is of order one and smaller, as expected for a well behaved perturbative expansion. 

2.2 High-beta simulations 

Calculation of higher order loop corrections becomes increasingly difficult due to the large 
number of diagrams and the complicated structure of vertices with a large number of gluons. 
Instead we can calculate radiative corrections by measuring the heavy quark two point function 
at large values of j3, extract the renormalization parameters and fit them to a polynomial in 
[14]. Nonperturbative contributions are suppressed by using twisted boundary conditions [ |T0| , p^ ]. 

'We only show the infrared finite part without the logarithmic divergence — 2/(3;r)logfl^A^. 
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Figure 2: One loop corrections to heavy quark renormalization parameters. Mean field corrections are 
included. The heavy quark mass is m = 2.8 and the stability parameter n = 2. 



fit 




v = 0.0 


v = 0.2 


v = 0.4 


v = 0.8 


unconstrained 


1 loop 


-2.6321(69) 


-2.5797(67) 


-2.4146(64) 


-1.812(16) 




2 loops 


-0.24(36) 


-0.33(36) 


-0.93(34) 


-0.01(67) 




3 loops 


-15.1±3.3 


-14.8±3.3 


-8.9±3.3 


-20.0 ±5.2 


constrained 


1 loop 


-2.6254 


-2.5732 


-2.4140 


-1.7284 




2 loops 


-0.55(17) 


-0.63(17) 


-0.96(15) 


3.06(59) 




3 loops 


-12.7±2.1 


-12.4±2.3 


-8.6±2.1 


0.3 ±5.9 



Table 1: Coefficients of the perturbative expansion of the zero point energy shift Eq for a simple mNRQCD 
action with 5H = 0. The calculation was carried out on a 18x6^ lattice and results are shown both for 
unconstrained and constrained one loop coefficients. Mean field corrections are not included. The heavy 
quark mass is m = 2.0 and the stability parameter n = 2. All results are preliminary. 

As the inclusion of fermionic vacuum polarisation effects makes the generation of configurations 
computationally expensive, the configurations used in our calculation are quenched. The missing 
fermionic corrections can be included relatively easily as they require the evaluation of a small 
number of diagrams. In addition, the quality of the polynomial fit can be improved if the one loop 
coefficient is constrained to the value calculated in diagrammatic perturbation theory. 

In Tab. |l] we show the results of a polynomial fit in a^. In the first case, the one loop coefficient 
was unconstrained, whereas in the second case it was constrained to the value from diagrammatic 
perturbation theory. For all but the highest frame velocity the one loop coefficients agree within 
statistical errors. 

2.3 Higher order kinetic terms 

We also investigate the breaking of rotational invariance by higher order kinetic terms in the 
NRQCD Lagrangian. Energy splittings between mesons with the same but different Y,j P^j have 
been calculated nonperturbatively in NRQCD. Although these splittings are small compared to the 
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Figure 3: One loop renormalization parameters of higher order kinetic terms in the NRQCD action. Mean 
field corrections are included. The stability parameter is n = 4. 



total kinetic energy, they are not compatible with zero within errors. They can be removed by 
tuning the coefficient of the rotationally noninvariant term A'"*' / (24m) to cs = 2.6, which has to be 
compared to the tree level value of c^^^^ = 1 . 

In Fig. ^ the one loop corrections in the perturbative calculation are shown, together with the 
corresponding corrections to the coefficient ci which multiplies the operator — (A'^') / (8m^). As 
can be seen from this plot, the smallness of the one loop corrections is not compatible with the 
nonperturbative calculation. Effects of higher order derivative operators in the action and the size 
of discretization errors in the nonperturbative simulations are currently being explored. 

3. Radiative corrections to heavy-light currents 

We calculate the one loop correction to the vector {V) and tensor (T) current. 

(2(^)^^ = ^y/^^F, e^^^^^ = -i^m^a'^^*^. (3.1) 

loTT^ 

The branching ratio of the decay B — > n£v can be calculated if the form factor of the vector current 
is known. The local contributions to the rare decays B K*y and B K^'*^U are given by the 
tensor and vector current respectively. 

The QCD field *F has to be expressed in terms of fields in the effective theory in the moving 
frame. At leading order in the 1/m expansion two operators Q\. contribute on the lattice for each 
Dirac structure P. From these a lattice operator Q\^^ can be constructed as 

Q^i = {l + a4^)QV + aJ'^d'\ (3.2) 
(r) 

The matching coefficients Cj_ are adjusted such that matrix elements in lattice mNRQCD and 
continuum QCD agree at one loop, {s\Q^^'i [^)qcd = (■^iGi^ !^)iat + ^{ccj, 1 />n). This implies 

c'^P{am,n/m) = 8Z^^\^l/m) - 8z'fl{am), c^^\am) = -5zJ2(am). (3.3) 
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frame velocity v frame velocity v 



Figure 4: One loop matching coefficients for the vector (left) and tensor (right) current. The Lorentz indices 
can be timelike (0), parallel (||) or perpendicular (IS) to the frame velocity. Mean field corrections are 
included, the massless light quark is discretized in the ASQtad action. The heavy quark mass is m 2.8 
and the stability parameter n = 2; the renormalization scale for the tensor current is jj. = m. 



the one loop correction to the continuum operator Q^^"^ whereas the one loop mixing 
matrix between the operators Q^f^ on the lattice is 5Z^^\am). 

Numerical results for the matching coefficients of the vector and tensor current are shown in 
Fig. ^ The light quark is discretized using the ASQtad action and its mass is set to zero in the 
matching calculation. Nonperturbative form factors of these currents have been calculated for a 
simple mNRQCD action [jl^]. This calculation is currently repeated with the full i^{l/m^,vf^^) 
mNRQCD action used in this work [19], including 1/m corrections to the currents at tree level. 



4. Conclusion and Outlook 

We calculated radiative corrections to leading order kinetic terms in the mNRQCD action on 
the lattice in diagrammatic perturbation theory both in an infinite volume and on a finite lattice 
with twisted boundary conditions. The one loop coefficients are compared to those from high-j3 
simulations and used to stabilise the polynomial fit in the strong coupling constant. We find good 
agreement between the two methods for a reasonable range of frame velocities. 

Radiative corrections to rotationally noninvariant higher order kinetic terms are calculated in 
NRQCD. Further work is necessary to reduce the mismatch of one loop results and nonperturbative 
calculations of the matching coefficients. Results for the one loop corrections to the leading order 
operators of the heavy-light vector and tensor currents in mNRQCD are presented. These will be 
combined with the nonperturbative calculation of the form factors which is currently being carried 
out with the full v^gj) action used in perturbative calculations. 

A two loop calculation of the fermionic contribution to the heavy quark self-energy is cur- 
rently being carried out and will complement the gluonic contributions from quenched high-j3 
simulations. Other terms in the action, most importantly the a ■ B/{2m) chromomagnetic term 
which affects the hyperfine splitting, need to be renormalized. Our approach for evaluating Feyn- 
man diagrams is flexible enough to extend it to other actions; so far we have repeated the current 
matching calculation with the HISQ light quark action but other combinations such as static heavy 
quarks decaying to domain wall fermions could be investigated as well. 
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